Search for: All records

ABSTRACT We introduce a new approach for analysing the intergalactic medium (IGM) damping wings imprinted on the proximity zones of quasars in the epoch of reionization (EoR). Whereas past work has typically forgone the additional constraining power afforded by the blue side continuum ($$\lambda \lesssim 1280\,$$ Å) and/or opted not to model the large correlated IGM transmission fluctuations in the proximity zone ($$\lambda \lesssim 1216\,$$ Å), we construct a generative probabilistic model for the entire spectrum accounting for all sources of error – the stochasticity induced by patchy reionization, the impact of the quasar’s ionizing radiation on the IGM, the unknown intrinsic spectrum of the quasar, and spectral noise. This principled Bayesian method allows us to marginalize out nuisance parameters associated with the quasar’s radiation and its unknown intrinsic spectrum to precisely measure the IGM neutral fraction, $$\langle x_{\rm H\,\small{I}}\rangle$$. A key element of our analysis is the use of dimensionality reduction (DR) to describe the intrinsic quasar spectrum via a small number of nuisance parameters. Using a large sample of 15 559 SDSS/BOSS quasars at $$z \gtrsim 2.15$$ we trained and quantified the performance of six distinct DR methods, and find that a six parameter principal component analysis model (five coefficients plus a normalization) performs best, with complex machine-learning approaches providing no advantage. By conducting statistical inference on 100 realistic mock EoR quasar spectra, we demonstrate the reliability of the credibility contours that we obtain on $$\langle x_{\rm H\,{\small{I}}}\rangle$$ and the quasar lifetime, $$t_{\rm Q}$$. The new method introduced here will transform IGM damping wings into a precision probe of reionization, on the same solid methodological and statistical footing as other precision cosmological measurements. 

ABSTRACT We present optical and near-infrared (NIR) spectroscopic observations for a sample of 45 quasars at $$6.50 < z \le 7.64$$ with absolute magnitudes at 1450 Å in the range $$-28.82 \le M_{1450} \le -24.13$$ and their composite spectrum. The median redshift and $$M_{1450}$$ of the quasars in the sample are $$z_{\rm {median}}=6.71$$ and $$M_{1450,\rm {median}} \simeq -26.1$$, respectively. The NIR spectra are taken with Echelle spectrographs, complemented with additional data from optical long slit instruments, and then reduced consistently using the open-source Python-based spectroscopic data reduction pipeline PypeIt. The median of the mean signal-to-noise ratios per 110 km s$$^{-1}$$ pixel in the J, H, and K band [median $$\langle \rm {SNR}_{\lambda } \rangle$$] is median $$\langle \rm {SNR}_{J} \rangle =9.7$$, median $$\langle \rm {SNR}_{H} \rangle =10.3$$, and median $$\langle \rm {SNR}_{K} \rangle =11.7$$; demonstrating the good data quality. This work presents the largest medium-/moderate-resolution sample of quasars at $z>6.5$ from ground-based instruments. Despite the diversity in instrumental set-ups and spectral quality, the data set is uniformly processed and well-characterized, making it ideally suited for several scientific goals, including the study of the quasar proximity zones and damping wings, the Ly $$\alpha$$ forest, the intergalactic medium’s metal content, as well as other properties such as the distribution of SMBH masses and Eddington ratios. Our composite spectrum is compared to others at both high and low z from the literature, showing differences in the strengths of many emission lines, probably due to differences in luminosity among the samples, but a consistent continuum slope, which proves that the same spectral features are preserved in quasars at different redshift ranges. 

ABSTRACT Lyman α emitters (LAEs) are excellent probes of the reionization process, as they must be surrounded by large ionized bubbles in order to be visible during the reionization era. Large ionized regions are thought to correspond to overdense regions and may be protoclusters, making them interesting test-beds for early massive structures. Close associations containing several LAEs are often assumed to mark overdense, ionized bubbles. Here, we develop the first framework to quantify the ionization and density fields of high-z galaxy associations. We explore the interplay between (i) the large-scale density of a survey field, (ii) Poisson noise due to the small number density of bright sources at high redshifts (z ∼ 7), and (iii) the effects of the ionized fraction on the observation of LAEs. We use Bayesian statistics, a simple model of reionization, and a Monte Carlo simulation to construct a more comprehensive method for calculating the large-scale density of LAE regions than previous works. We find that Poisson noise has a strong effect on the inferred density of a region and show how the ionized fraction can be inferred. We then apply our framework to the strongest association yet identified: Hu et al. found 14 LAEs in a volume of ∼50 000 cMpc3 inside the COSMOS field at z ∼ 7. We show that this is most likely a 2.4σ overdensity inside of an ionized or nearly ionized bubble. We also show that this LAE association implies that the global ionized fraction is $$\bar{Q} = 0.59^{+0.10}_{-0.11}$$, within the context of a simple reionization model. 

ABSTRACT We present the first observational measurements of the Lyman-α (Ly α) forest flux autocorrelation functions in ten redshift bins from 5.1 ≤ z ≤ 6.0. We use a sample of 35 quasar sightlines at z > 5.7 from the extended XQR-30 data set; these data have signal-to-noise ratios of >20 per spectral pixel. We carefully account for systematic errors in continuum reconstruction, instrumentation, and contamination by damped Ly α systems. With these measurements, we introduce software tools to generate autocorrelation function measurements from any simulation. Our measurements of the smallest bin of the autocorrelation function increase with redshift when normalizing by the mean flux, 〈F〉. This increase may come from decreasing 〈F〉 or increasing mean free path of hydrogen-ionizing photons, λmfp. Recent work has shown that the autocorrelation function from simulations at z > 5 is sensitive to λmfp, a quantity that contains vital information on the ending of reionization. For an initial comparison, we show our autocorrelation measurements with simulation models for recently measured λmfp values and find good agreements. Further work in modelling and understanding the covariance matrices of the data is necessary to get robust measurements of λmfp from this data. 

Abstract The [Cii] 158μm emission line and the underlying far-infrared (FIR) dust continuum are important tracers for studying star formation and kinematic properties of early galaxies. We present a survey of the [Cii] emission lines and FIR continua of 31 luminous quasars atz> 6.5 using the Atacama Large Millimeter Array (ALMA) and the NOrthern Extended Millimeter Array at sub-arcsec resolution. This survey more than doubles the number of quasars with [Cii] and FIR observations at these redshifts and enables statistical studies of quasar host galaxies deep into the epoch of reionization. We detect [Cii] emission in 27 quasar hosts with a luminosity range ofL_[CII]= (0.3–5.5) × 10⁹L_⊙and detect the FIR continuum of 28 quasar hosts with a luminosity range ofL_FIR= (0.5–13.0) × 10¹²L_⊙. BothL_[CII]andL_FIRare correlated (ρ≃ 0.4) with the quasar bolometric luminosity, albeit with substantial scatter. The quasar hosts detected by ALMA are clearly resolved with a median diameter of ∼5 kpc. About 40% of the quasar host galaxies show a velocity gradient in [Cii] emission, while the rest show either dispersion-dominated or disturbed kinematics. Basic estimates of the dynamical masses of the rotation-dominated host galaxies yieldM_dyn= (0.1–7.5) × 10¹¹M_⊙. Considering our findings alongside those of literature studies, we found that the ratio betweenM_BHandM_dynis about 10 times higher than that of localM_BH–M_dynrelation on average but with substantial scatter (the ratio difference ranging from ∼0.6 to 60) and large uncertainties. 

ABSTRACT Recent quasar absorption line observations suggest that reionization may end as late as $$z \approx 5.3$$. As a means to search for large neutral hydrogen islands at $$z\ \lt\ 6$$, we revisit long dark gaps in the Ly $$\beta$$ forest in Very Large Telescope/X-Shooter and Keck/Echellette Spectrograph and Imager quasar spectra. We stack the Ly $$\alpha$$ forest corresponding to both edges of these Ly $$\beta$$ dark gaps and identify a damping wing-like extended absorption profile. The average redshift of the stacked forest is $z=5.8$. By comparing these observations with reionization simulations, we infer that such a damping wing-like feature can be naturally explained if these gaps are at least partially created by neutral islands. Conversely, simulated dark gaps lacking neutral hydrogen struggle to replicate the observed damping wing features. Furthermore, this damping wing-like profile implies that the volume-averaged neutral hydrogen fraction must be $$\langle x_{\rm H\,{\small {I}}} \rangle \ge 6.1 \pm 3.9~{{\ \rm per\ cent}}$$ at $z = 5.8$. Our results offer robust evidence that reionization extends below $z=6$.